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Изобретение относится к измерительной 
технике и может найти применение в 
навигационных системах. Технический результат - 
повышение точности измерений. Для достижения 
данного результата измеряют геодезические 
координаты размещенного на носителе 
навигационного комплекса и с помощью 
размещенного там же устройства определения 
положения объекта - вектор положения объекта 
относительно устройства. Определяют параметры 
вектора положения объекта относительно 
устройства и вектора положения навигационного 
комплекса. При этом дополнительно измеряют 


вектор положения устройства относительно 


навигационного комплекса и угловую ориентацию 
носителя относительно земли. 8 з.п. ф-лы, 1 ил. 
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(54) MODE OF DEFINITION OF THE LAND-SURVEYING COORDINATES OF AN OBJECT 


(57) Abstract: 

FIELD: the invention refers to measuring 
technique and may find application in 
navigational systems. 

SUBSTANCE: for achieving the given result land- 
surveying coordinates of the navigational complex 
placed on a carrier are measured and with the aid 
of an arrangement there placed of definition of 
the object relatively to the arrangement. The 
parameters of the vector of the position of the 
object relatively to the arrangement and_ the 
vector of the position of the navigational 
complex are defined. At that the vector of the 
position of the arrangement relatively to the 


navigational complex and the angular 


orientation 


of the carrier relatively to the Earth are defined. 
EFFECT: increases accuracy of measurements. 
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Изобретение относится к геодезии, в частности к способам определения геодезических 
координат объекта. Объектами могут служить участки земной поверхности, строения, 
воздушные и космические летательные аппараты (JIA), морские и речные суда ит.п. 
Изобретение может быть использовано в системах воздушной, морской и космической 
навигации, для целеуказания, а также при наземной послеполетной обработке 
зарегистрированных данных средств бортовых измерений летательного аппарата для 
анализа погрешностей определения геодезических координат объекта с известными 
геодезическими координатами бортовыми средствами, в частности радиолокатором, 
лазерным дальномером. 

Для определения геодезических координат объекта на носителе (в качестве которого 
может быть использован летательный аппарат, корабль, наземная стационарная или 
подвижная установка) размещают навигационный комплекс и устройство измерения 
относительного положения D объекта (т.е. положения объекта относительно этого 
устройства). При этом геодезические координаты навигационного комплекса (долгота La, 
широта Ba и высота Ha) либо определены заранее путем выполнения его геодезической 
привязки (наземные стационарные радиосистемы ближней или дальней навигации), либо 
измеряются на борту подвижного носителя с помощью, например, глобальной 
навигационной спутниковой системы (ГНСС) [ГОСТ Р 51794-2001. Государственный 
стандарт Российской Федерации. Аппаратура радионавигационная глобальной 
навигационной спутниковой системы и глобальной системы позиционирования. Системы 
координат. Методы преобразований координат определяемых точек. - М.: Госстандарт 
России, с.1...2]. В дальнейшем рассмотрим, как общий, случай размещения 
навигационного комплекса и устройства измерения положения объекта на борту 
подвижного носителя. 

Определение геодезических координат: широты В, долготы L и высоты Н объекта по 
измеренным геодезическим координатам: широте Ba, долготе La и высоте Нд 
навигационного комплекса и измеренному вектору относительной дальности р до 
объекта, называется прямой пространственной геодезической задачей и представлено, 
например, в [Морозов В.П. Курс сфероидической геодезии. - М.: Недра, 1979, с.196-197]. 

Определим используемую ниже терминологию. Вектор положения какого-либо объекта в 
референцной геодезической системе координат (СК), связанной с общеземным 
эллипсоидом (O39), задают геодезическими координатами: широтой В, долготой L n 
высотой Н [ГОСТ Р 51794-2001, с.1...2]. Размер и форму O39 задают двумя параметрами: 
большой полуосью а и эксцентриситетом е [там же, с.1, 3]. Геоцентрической широтой Ф 
объекта называют угол между экваториальной плоскостью O39 и прямой, соединяющей 
центр ОЗЭ с объектом [Морозов В.П. Курс сфероидической геодезии. - М.: Недра, 1979, 
с.189]. 

Под угловой ориентацией подвижного носителя навигационного комплекса (например, 
ЛА) относительно земли понимают значения углов Эйлера, в частности тангажа 3, крена y 
и рыскания y [ГОСТ 20058-80. Государственный стандарт Союза ССР. Динамика 
летательных аппаратов в атмосфере. Термины, определения и обозначения. - М.: 
Государственный комитет СССР по стандартам, 1981. с.9], либо матрицу направляющих 
косинусов [там же, с.43] между осями OX, OY, OZ связанной СК и осями ОХ,, ОҮ,, ОД, 
нормальной СК [там же, с.4]. Начало О нормальной СК OX, Yg 2, совпадает с началом 
связанной СК подвижного носителя, ось OY, направлена по нормали к поверхности O39, 
ось ОХ. направлена в сторону оси вращения ОЗЭ. Угол курса y, отличается от угла 
рыскания y противоположным направлением положительного отсчета. Как видим, 
нормальная СК, принятая в динамике ЛА в атмосфере, совпадает с топоцентрической 
декартовой прямоугольной СК Opgcy, принятой в сфероидической геодезии. Вектор 
относительной дальности р до объекта в топоцентрической СК задают полярными 
координатами - наклонной относительной дальностью D, геодезическим азимутом Аб и 
зенитным расстоянием Zo [Морозов В.П. Курс сфероидической геодезии. - M.: Недра, 
1979, с.189...190]. 
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Введем локальную декартову геоцентрическую прямоугольную СК O,X-Y,Z,, связанную 
с положением подвижного носителя навигационного комплекса. Начало О, ее поместим в 
центре ОЗЭ. Плоскость O,X,Y, расположим в меридианной плоскости, содержащей начало 
О нормальной СК ЛА, ось ОХ; направим вдоль пересечения меридианной плоскости с 
экваториальной плоскостью O39 от точки O, к проекции точки О на это пересечение, 
ось ОгҮ; направим в меридианной плоскости на северный полюс ОЗЭ, ось O,Z, дополняет 


СК до правой и лежит в экваториальной плоскости ОЗЭ. 


Прототипом изобретения является способ определения геодезических координат 
объекта, представленный в [Морозов В.П. Курс сфероидической геодезии. - М.: Недра, 
1979, с.196...197], в котором измеряют геодезические координаты: долготу La, широту Вди 
высоту Нд размещенного на носителе навигационного комплекса и с помощью 
размещенного на упомянутом носителе устройства определения положения объекта - 
вектор положения р =(0 Ау Zo) объекта относительно упомянутого устройства, по 
измеренным данным вычисляют в декартовой геоцентрической системе координат 
компоненты вектора положения Ё, = (Rax, Ray Raz) навигационного комплекса и 


вектора положения р = (р, р Be ) объекта относительно упомянутого устройства и 
г г ёг 


определяют в этой системе координат вектор положения È = (R R 


Хг 1 


R, объекта 
p 


суммированием вычисленных упомянутых компонент векторов, по вектору положения 
объекта определяют геодезические координаты объекта: долготу L, широту В и высоту Н. 


Компоненты Ray Raye Raz, 


декартовой геоцентрической СК вычисляют по формулам 


Е ee 
J1- е sin’ В, 

Rax = (№, +H, )соѕВ, , (1) 

Кар = la ` eĉ Na + н, [зіп Bay 

70, 


где Na - вспомогательная переменная, а компоненты D gnr |0 


р объекта относительно устройства вычисляют по формулам 


р 
p7 “Zp 


вектора положения Ra навигационного комплекса в 


вектора положения 


Dey = DSim COS Aes 
Dy, = Dcosdé,, 
De, = 05111, SinA,, 
Dy = Dy, Ba Baot Dy, COSB; 
Dy = Dy, cosB, + Dy ВІВ у 
Da, = =. 
roe Dyg Dy s De, - компоненты вектора положения р объекта в нормальной CK. 


В известном способе для определения геодезической долготы L объекта сначала 


вычисляют приращение | геодезической долготы объекта 


R 
| = arctg —^ 2 
агі (2) 


х, 


по компонентам Ry, и А, вектора положения В, объекта по осям геоцентрической 


системы координат в экваториальной ПЛОСКОСТИ общеземного эллипсоида с учетом знаков 


этих компонентов, а затем саму геодезическую долготу объекта 
= А+. 


Геодезическую широту В объекта определяют из уравнения связи тангенса ї 
геодезической широты с компонентами Ry; Ry , А, вектора положения В объекта 
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ee ee E С. (3) 
vk +t? 


в котором введены вспомогательные переменные и константы: 


Е = typ (4) 
2 2 
Bip = Вы ов: 
Ry 
t, = = 5 
p 9 
= 1 
bes 
аке? 
ро 
R h 


Как видно из правой части выражения (5), величина to представляет значение тангенса 
геоцентрической широты объекта. 
Для вычисления геодезической широты объекта применяют итерационную процедуру 


= (ais Oi ge css) (6) 
qk + ti 
до достижения заданной точности значения тангенса t геодезической широты В. Из (6) 
видно, что геодезическую широту объекта определяют многократным уточнением 
приближенного значения ее тангенса непосредственно из уравнения связи (3). После чего 
из выражения (4) определяют саму геодезическую широту объекта 
В = arctgt,,;} (7) 
Геодезическую высоту Н объекта вычисляют по формуле 


Ешан а = (3) 


h 5 1+1 
yk T tang 


Недостатком этого способа определения геодезических координат L, В, H объекта 
являются неприемлемые для некоторых целей погрешности ёі, 5B, Н их определения. 
Эти погрешности вызваны следующими обстоятельствами. 

Во-первых, измеряемые геодезические координаты навигационного комплекса носителя 
являются координатами его определенной части. Так, при использовании на носителе 
аппаратуры [HCC измеряются координаты La, Ba n Ha антенны этой аппаратуры. Вектор 
положения D объекта определяет вектор положения объекта относительно размещенного 
на упомянутом носителе устройства определения положения объекта. Антенна аппаратуры 
ГНСС упомянутого навигационного комплекса и упомянутое устройство пространственно 
разнесены относительно друг друга, причем величина этого разноса может изменяться от 
нескольких метров на ЛА до десятков и сотен метров на корабле и определяет величину 
методической погрешности измерения вектора относительного положения объекта в 
известном способе. 

Обозначим в декартовой геоцентрической СК через Е, вектор положения 
навигационного комплекса, через г - вектор положения упомянутого устройства 
относительно упомянутого навигационного комплекса, через Г - вектор положения 
объекта относительно упомянутого устройства, тогда вектор положения Е объекта 
определяется следующей суммой (фиг.1): 

В = В, +Е+6. (9) 

В известном способе определения геодезических координат объекта вектор его 
положения R, определяют согласно выражениям (1) другой суммой 
Ви = К, +5 (10) 

не содержащей вектора положения г объекта относительно упомянутого устройства. 
Это приводит в известном способе к погрешности 58, вычисления геоцентрических 


Е 


itl 
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декартовых координат объекта, определяемой из выражений (9) и (10) 


— 


SR, = Ry - Е = - (11) 

Максимальные погрешности измерения дальности современными лазерными 
дальномерами не превышают 1 м. Следовательно, указанные выше значения (от 
нескольких до десятков и сотен метров) погрешностей 58, вычисления геоцентрических 
декартовых координат объекта известным способом могут оказаться недопустимо 
большими. Это приводит к соответствующим погрешностям 6, 6B, SH определения 
геодезических координат объекта. Поэтому известный способ может приводить к 
определению геодезических координат объекта с недопустимо большой методической 
погрешностью. 

Во-вторых, снижение погрешности $В вычисления геодезической широты В объекта при 
итерационной процедуре (6) проходит медленно и с неизвестным заранее количеством 
итераций (для достижения заданной точности). В вычислительных системах реального 
времени неопределенное количество итераций является существенным недостатком. 

В-третьих, в известном способе при определении приращения | (2) геодезической 
долготы объекта отношение А, неограниченно возрастает по модулю при значении 


Вх, 


Вх, = 0 , Которое может выполняться при любых значениях геодезической широты В 


объекта. Поэтому вычисление значений тангенсов как геодезической долготы | (2), так и 
геодезической широты В (7) объекта приводит к резкому возрастанию погрешностей их 
расчета при значениях углов, близких к Г д. 


Е. 
Действительно, из выражения (7) следует, что погрешности $В измерения геодезической 
широты объекта приводят к следующим погрешностям вычисления значения ее тангенса 
ôE 


cos E | 
Правая часть данного выражения неограниченно возрастает при значениях 
геодезической широты объекта, близких к т. Правая часть выражения (2) для расчета 
0 
приращения | геодезической долготы объекта содержит аналогичную особенность при 
значении к. = 0 составляющей Ry, вектора положения объекта, т.е. при значениях 


öt = 


приращения | геодезической долготы объекта, близкихк т. 
2 
В-четвертых, указанное резкое возрастание погрешности расчета тангенса 
геодезической широты В (4) при значениях углов, близких к + п, согласно выражению (8) 
5 
приводит к неопределенному результату расчета геодезической высоты Н объекта, 


поскольку при этом величина Ry стремится к нулю, множитель l es K 


бесконечности, а их произведение может дать произвольный непредсказуемый результат. 

Задачей изобретения является повышение точности определения геодезических 
координат объекта. 

Задача решается с помощью способа определения геодезических координат объекта, в 
котором измеряют геодезические координаты размещенного на носителе навигационного 
комплекса и с помощью размещенного на упомянутом носителе устройства определения 
положения объекта - вектор положения объекта относительно упомянутого устройства, по 
измеренным данным вычисляют в декартовой геоцентрической системе координат 
компоненты вектора положения объекта относительно упомянутого устройства и вектора 
положения навигационного комплекса и определяют в этой системе координат вектор 
положения объекта суммированием вычисленных упомянутых компонент, по вектору 
положения объекта определяют геодезические координаты объекта, отличающегося тем, 
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что дополнительно измеряют вектор положения упомянутого устройства относительно 
упомянутого навигационного комплекса и угловую ориентацию носителя относительно 
земли, при упомянутом вычислении компонент в декартовой геоцентрической системе 
координат дополнительно вычисляют компоненты упомянутого вектора положения 
упомянутого устройства относительно упомянутого навигационного комплекса и 
дополнительно их учитывают при упомянутом суммировании. 

С целью снижения погрешности определения геодезических широты и высоты объекта 
геодезическую широту объекта определяют уточнением как минимум один раз 
приближенного значения ее тангенса { путем вычитания корректирующей поправки, равной 
отношению значения функции 
РЕ 


vk + t? 


где to - тангенс геоцентрической широты объекта, 


F(t) = t - - 


k 1 
— о = 
1 – е“ 
е - эксцентриситет общеземного эллипсоида, 
aske? 
oS 
А 


h 

а - большая полуось общеземного эллипсоида, 

В; - модуль проекции вектора положения объекта на экваториальную плоскость 
общеземного эллипсоида, к производной этой функции по t при значении t, равном 
предшествующему приближенному его значению. 

С целью снижения погрешности определения геодезических широт начальное 
приближенное значение тангенса t геодезической широты определяют по тангенсу to 
геоцентрической широты объекта по формуле 


pty 


yk + to | 


С целью снижения погрешности определения геодезических широты и высоты объекта 
геодезическую широту В объекта определяют уточнением как минимум один раз 
приближенного ее значения путем вычитания корректирующей поправки, равной 
отношению значения функции 


sinBcosB 


Vi-e?sin?B | 

где  - геоцентрическая широта объекта, 

_ ae’ 
Я a 

R - модуль вектора положения объекта, 

к производной этой функции по В при значении В, равном предшествующему 
приближенному ее значению. 

С целью снижения погрешности определения геодезических широты и высоты объекта 
начальное приближенное значение геодезической широты В определяют по 
геоцентрической широте Ф объекта по формуле 


sin ФсоѕФ 


\1-е? sin? Ф 

С целью снижения погрешности определения геодезических широты и высоты объекта 
геоцентрическую широту Ф объекта определяют по величинам отношений проекции 
вектора положения объекта на ось вращения общеземного эллипсоида и модуля проекции 
этого вектора на экваториальную плоскость общеземного эллипсоида к модулю этого 
вектора. 

C целью снижения погрешности определения геодезической высоты объекта 


t=t, + 


F(B)=B- Ф агсѕіп| q 


В = Ф +агсѕіп| q 
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геодезическую высоту Н объекта определяют по формуле 
Н = В, созВ + В, 5118 - ayl - ef sin’ В, 
roe Ву. - компонент вектора положения объекта в декартовой геоцентрической системе 


координат, параллельный оси вращения общеземного эллипсоида. 

С целью снижения погрешности определения геодезической долготы объекта 
геодезическую долготу объекта определяют по величинам отношений тех компонентов 
вектора положения объекта в геоцентрической системы координат, которые соответствуют 
осям, лежащим в экваториальной плоскости общеземного эллипсоида, к модулю проекции 
данного вектора на упомянутую плоскость. 

Предлагаемое изобретение позволяет определять геодезические координаты объекта с 
более высокой точностью. 

Применение предлагаемого способа иллюстрируется чертежом, на котором в 
декартовой геоцентрической системе координат показаны векторы: 

В, - положения навигационного комплекса; 

г - положения устройства определения положения объекта относительно упомянутого 
навигационного комплекса; 

D - положения объекта относительно устройства определения положения объекта; 

в, - положения объекта в известном способе; 

В - положения объекта в предлагаемом способе; 

$, - погрешности определения положения объекта известным способом. 

Предлагаемый способ определения геодезических координат объекта заключается в 
следующем. 

1. Измеряют геодезические координаты La, Ba, Ha размещенного на носителе 
навигационного комплекса и с помощью размещенного на упомянутом носителе устройства 
определения положения объекта - вектор положения р = (р А; Za) объекта 


относительного упомянутого устройства. 

2. Дополнительно измеряют вектор положения É = (к; ry г, упомянутого 
устройства относительно упомянутого навигационного комплекса и угловую 
ориентацию 5, ү, ли, носителя относительно земли. 

Измерение вектора положения ғ упомянутого устройства относительно упомянутого 
навигационного комплекса и угловой ориентации 5, у, ли, носителя относительно земли 
является существенным признаком изобретения, поскольку необходимо для исключения 
методической погрешности (11) определения геодезических координат объекта. 

3. По измеренным данным вычисляют в декартовой геоцентрической системе координат 


компоненты Raye Raye Каз, вектора положения Ra навигационного комплекса 
Ray, = (Ny + Hy) cos By; 


р = И - е? м, + H,|sin Bar 


где а - вспомогательная переменная, 


А 
| ree 
l1- esin B, 
и в этой же системе координат вычисляют компоненты Dy , Dy, Dz, Вектора 


положения D объекта относительно упомянутого устройства 


Dy, = Dx, ВВ = Dy, COSB,, 
Dy, = Ds. eos Ber Dy, SL Big 
Da, = Da? 
nA) = 05111, cosA, ‚В, = 25054, "Da = 051142, sinA, 


компоненты упомянутого вектора положения р объекта в нормальной СК. 
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4. Дополнительно В декартовой геоцентрической системе координат вычисляют 


компоненты ү, ‚Е, 6 вектора положения Е упомянутого устройства относительно 
T г г 
упомянутого навигационного комплекса 
г, = -Pry sinB, + ry, соБѕВ,, 
гү = ty CosB, + ry, SIN Bar 
г, = tay 
roe 
ry, = Ty СОЗ, cosg — г,(5іпц, Siny + cosw, 51135051) = r(siny, cosy — cosy, 51135117), 
ty, = 1; 5113 + Е, С053 соѕу – г, cos$siny, 


lz 
положения у В нормальной системе координат. 
5. Определяют в декартовой геоцентрической системе координат вектор положения 


В = (Ах, Ку Re) объекта суммированием вычисленных упомянутых компонент 


векторов Ё, = (Ray Ray, Raz,) ИВ (ye D 
компоненты вектора г = (ry г, г, ) положения устройства относительно 


= г, siny,cos$+r, (cosy, siny—siny, sin $cosy)+r (cosy, cosy + siny, sin9siny)- компоненты вектора 


Dz, ) ‚ дополнительно учитывая 


навигационного комплекса: 
Ry =R,y + Dy + Ky 

Be St Bae ВЕ ЫЕ, 
Be Во Da ЯА 


6. По вектору положения Е = (в, R, R объекта определяют геодезические 
T T 


z) 
координаты объекта: долготу L, широту B и высоту H. 

С целью снижения погрешности определения геодезической широты В;.; объекта при 
определении геодезических координат объекта по вектору положения объекта 
геодезическую широту находят путем приближений ее значений. Предлагаются 2 варианта. 

Вариант 1. Геодезическую широту В объекта определяют уточнением как минимум один 
раз приближенного значения ее тангенса t путем вычитания корректирующей поправки At, 
т.е. 


АЕ Е (12) 
Поправка At равна отношению значения функции 
t 
Ре. (13) 


vk +t? 


к скорости изменения этой функции (т.е. к производной этой функции по t), при 
значении t, равном предшествующему приближенному его значению +, т.е. 


м (м) 


[а] 


Функция F(t) представляет собой левую часть уравнения (3), связывающего тангенс t 
геодезической широты с декартовыми геоцентрическими координатами объекта (4)...(5). 

Уточнение (12) приближенного значения тангенса ti путем введения корректирующей 
поправки (14) основано на численном решении нелинейного уравнения методом Ньютона 
[Самарский A.A., Гулин А.В. Численные методы: Учебное пособие - M.: Наука. Гл. ред. 
физ.-мат. лит., 1989, с.193]. Известно [там же, с.200], что сходимость метода Ньютона 
существенно более быстрая, чем метод простой итерации (6), использованный в известном 
способе. 

Начальное приближенное значение тангенса t; геодезической широты определяют по 
тангенсу tg геоцентрической широты объекта из уравнения (3) по формуле 
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t 
pobo 


i 0 м z 
Jeet 
Данное определение начального приближенного значения тангенса геодезической 
широты обеспечивает более высокую точность определения геодезической широты 
объекта вследствие использования уравнения связи (3) геодезической широты объекта с 
компонентами вектора положения объекта в декартовой геоцентрической СК. 
Вариант 2. Геодезическую широту В определяют уточнением как минимум один раз ее 
приближенного значения путем вычитания корректирующей поправки AB 
В..1 = В, - AB (15) 
Поправка AB пропорциональна отношению значения функции 
sinBcosB \ 6) 
м1 е? sinB 
к производной этой функции по В (скорости изменения dF (в 


ав 
равном предшествующему приближенному ее значению, т.е. 


F(B)=B- ®-aresin| q 


) ) при значении B, 


Agee). Gz) 


Выражение (16) получено из уравнения связи геодезической широты с декартовыми 
геоцентрическими координатами объекта при условии 
F(B)=0. (18) 

Операция (15) уточнения приближенного значения геодезической широты В; путем 
введения корректирующей поправки AB (17) также представляет собой численное решение 
нелинейного уравнения (18) методом Ньютона. Уравнение (18) имеет следующий вид: 


sinBcosB 
В-Ф- arcsin| а |= 0, (19) 
41 е? ѕіп В 
где 
2 a 
=e —, 
i В 
В = В, +В, "+В, ", 
R 
їг 
Ве = 2 20 
2-3 20 
Ry 
ey eo Zi 
pa (21) 
arcsins при Cp = Е 
7 a 
Bo ИЯ 
Ф - 2 (22) 
= — — āarcsince ODPM зое 
Bq 2 
- + arcsinc при 3. < - 
2 $ Eg 7 Ia 
Здесь R - модуль вектора положения объекта в декартовой геоцентрической системе 
координат. 


Уравнение (19) получено преобразованием уравнения (3) с целью снижения 
погрешности определения геодезической широты В объекта (t стремится по модулю к 
бесконечности при стремлении В к 2 ). Это уравнение не содержит особенностей BO 

ВЕ 
всем диапазоне возможных значений геодезической широты. 
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Использование функции (16) для уточнения приближенного значения геодезической 
широты В; путем введения корректирующей поправки AB (17) обеспечивает повышение 
точности определения геодезической широты объекта. 

Начальное приближенное значение В; геодезической широты определяют по 
геоцентрической широте Ф объекта из уравнения (19) по формуле 


sin ФсоѕФ (23) 


41- e° sin? Ф 

Данное определение начального приближенного значения геодезической широты 
обеспечивает более высокую точность определения геодезической широты объекта 
вследствие использования уравнения связи (19) геодезической широты объекта с 
компонентами вектора положения объекта в декартовой геоцентрической системе 
координат. 

В предлагаемом способе геоцентрическую широту Ф объекта определяют согласно 
(20)...(22) по величинам отношений проекции вектора положения объекта на ось вращения 
общеземного эллипсоида (компонента Ry ) и модуля Rp проекции этого вектора Ha 


В, = ®+arcsin| q 


экваториальную плоскость общеземного эллипсоида к модулю К этого вектора. 
Известное более простое определение геоцентрической широты объекта Ф, например 
только по отношению проекции вектора положения объекта на вертикальную ось R, 
| 


декартовой геоцентрической СК к модулю К вектора положения объекта (20) 
®=arcsins,, (24) 

приводит к неограниченному возрастанию погрешностей §P ее вычисления вблизи 
значений т. Действительно, из выражения (24) следует для погрешности 


759 
геоцентрической широты 
05 
t =- В, (25) 


fl — зБ" 


В свою очередь согласно выражению (20) справедливо 

без = oR oe (26) 
R в? 

Из выражений (25) и (26) следует, что малейшие погрешности определения вектора 
положения Е объекта приводят к неограниченному возрастанию погрешности определения 
геоцентрической широты Ф объекта, а следовательно, согласно (23) и его геодезической 
широты B;. Поэтому предлагаемое определение геоцентрической широты Ф объекта 
обеспечивает повышение точности определения геодезической широты объекта. 

С целью снижения погрешности геодезическую высоту Н объекта определяют по 


формуле 
Н = В, cosB + Ry sinB = ayl - ef sin’ В, (27) 
где Ry - компонент вектора положения объекта в декартовой геоцентрической системе 


координат по оси вращения общеземного эллипсоида [ГОСТ P 51794-2001. 
Государственный стандарт Российской Федерации. Аппаратура радионавигационная 
глобальной навигационной спутниковой системы и глобальной системы позиционирования. 
Системы координат. Методы преобразований координат определяемых точек. - М.: 
Госстандарт России, с.3...4]. 

Определение геодезической высоты Н по (27) обеспечивает повышение точности 
определения геодезической высоты объекта при любых возможных значениях 
геодезической широты В в отличие от известного способа. 

С целью снижения погрешности геодезическую долготу L объекта определяют 


величинами отношений тех компонентов Rg y Ва, вектора положения объекта в 
г 
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геоцентрической системы координат, которые соответствуют осям, лежащим в 
экваториальной плоскости общеземного эллипсоида, к модулю проекции Rp данного 
вектора на упомянутую плоскость. Эти 


R 
а (28) 
Ry 
Ry 
= —, 29 
=i R, (29) 
aresin|s| при [$ < L 
42 
1, = 1 Ц 


z- агоэт при |5} 


42 


l= А (30) 


27 — Ш при c; 20 и 30 

L=La +l. 

Предлагаемый способ (28)...(30) определения приращения | геодезической долготы 
объекта содержит особенность типа 0/0 только при значении R,=0 (28, 29), которое 
соответствует нахождению объекта над одним из полюсов ОЗЭ. Но в полюсах значение 
геодезической долготы изначально не определено, поэтому при значении R,=0 для 
приращения геодезической долготы объекта можно положить |=0. 

Определение приращения | геодезической долготы объекта по величинам отношений 
компонентов в, , Кь вектора положения объекта по осям геоцентрической системы 


координат в экваториальной плоскости общеземного эллипсоида к модулю К проекции 
данного вектора на упомянутую плоскость обеспечивает повышение точности определения 
геодезической широты объекта. 

Предлагаемое изобретение позволяет определять геодезические координаты объекта с 
более высокой точностью, поскольку исключает методические погрешности определения 
геодезических координат объекта известным способом. 


Формула изобретения 

1. Способ определения геодезических координат объекта, при котором измеряют 
геодезические координаты размещенного на носителе навигационного комплекса и с 
помощью размещенного на упомянутом носителе устройства определения положения 
объекта - вектор положения объекта относительно упомянутого устройства, по измеренным 
данным вычисляют в декартовой геоцентрической системе координат компоненты вектора 
положения объекта относительно упомянутого устройства и вектора положения 
навигационного комплекса и определяют в этой системе координат вектор положения 
объекта суммированием вычисленных упомянутых компонент, по вектору положения 
объекта определяют геодезические координаты объекта, отличающийся тем, что 
дополнительно измеряют вектор положения упомянутого устройства относительно 
упомянутого навигационного комплекса и угловую ориентацию носителя относительно 
земли, при упомянутом вычислении компонент в декартовой геоцентрической системе 
координат дополнительно вычисляют компоненты упомянутого вектора положения 
упомянутого устройства относительно упомянутого навигационного комплекса и 
дополнительно их учитывают при упомянутом суммировании. 

2. Способ по п.1, отличающийся тем, что геодезическую широту объекта определяют 
уточнением как минимум один раз приближенного значения ее тангенса t путем вычитания 
корректирующей поправки, равной отношению значения функции 
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pt 
yk + t ' 


где to - тангенс геоцентрической широты объекта, 


Е) = t-t, - 


Е 1 
= в. 
е - эксцентриситет общеземного эллипсоида, 
aykef 
ee 
Ry 


а - большая полуось общеземного эллипсоида, 

Rp - модуль проекции вектора положения объекта на экваториальную плоскость 
общеземного эллипсоида, 

к производной этой функции по t при значении t, равном предшествующему 
приближенному его значению. 

3. Способ по п.2, отличающийся тем, что начальное приближенное значение тангенса { 
геодезической широты определяют по тангенсу t, геоцентрической широты объекта по 


формуле 


4. Способ по п.1, отличающийся тем, что геодезическую широту В объекта определяют 
уточнением как минимум один раз ее приближенного значения путем вычитания 
корректирующей поправки, равной отношению значения функции 


sinBcosB 


Ле? эт? В | 

где Ф - геоцентрическая широта объекта, 

aef 
= 
R 

R - модуль вектора положения объекта, 

к производной этой функции по В при значении В, равном предшествующему 
приближенному ее значению. 

5. Способ по п.4, отличающийся тем, что начальное приближенное значение 
геодезической широты В определяют по геоцентрической широте Ф объекта по формуле 


sinPcas? 


41-е? sin’ > 

6. Способ по п.5, отличающийся тем, что геоцентрическую широту Ф объекта 
определяют по величинам отношений проекции вектора положения объекта на ось 
вращения общеземного эллипсоида и модуля проекции этого вектора на экваториальную 
плоскость общеземного эллипсоида к модулю этого вектора. 

7. Способ по любому из пп.1-6, отличающийся тем, что геодезическую высоту Н объекта 
определяют по формуле 


H = R, Оов В Ry SinB- ayl = e sin’ В, 


где Ry, - компонент вектора положения объекта в декартовой геоцентрической системе 


ЕВ) =В-Ф- агсѕіп q 


B= Ф +агсѕіп! q 


координат, параллельный оси вращения общеземного эллипсоида. 

8. Способ по любому из пп.1-6, отличающийся тем, что геодезическую долготу объекта 
определяют по величинам отношений тех компонентов вектора положения объекта в 
геоцентрической системе координат, которые соответствуют осям, лежащим в 
экваториальной плоскости общеземного эллипсоида, к модулю проекции данного вектора 
на упомянутую плоскость. 

9. Способ по n.7, отличающийся тем, что геодезическую долготу объекта определяют по 
величинам отношений тех компонентов вектора положения объекта в геоцентрической 
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системе координат, которые соответствуют осям, лежащим в экваториальной ПЛОСКОСТИ 
общеземного эллипсоида, к модулю проекции данного вектора на упомянутую плоскость. 
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